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ABSTRACT 

X 

\ Chandra observations of solar-like pre-main sequence (PMS) stars in the 

Orion Flanking Fields (age ~1 Myr) and NGC 2264 (~3 Myr) are compared 
with the results of the COUP survey of similar objects in the ONC (~0.5 Myr). 
The correlations between log Lx and mass found for PMS stars on convective 
tracks in these clusters are consistent with the relationships found for the ONC, 
indicating little change in the median values of either log Lx or log L^/ Lboiduring 
the first ~3-5 Myr of evolution down convective tracks. The fraction of stars with 
extreme values of Lx, more than 10 times higher than the average for a given 
Lboi or with log Lx/Lboi greater than the canonical saturation value of —2.9, is 
however larger by a factor of two in the younger ONC when compared with the 
Orion FF and NGC 2264. 

PMS stars in NGC 2264 on radiative tracks have Lx/Lboi values that are 
systematically lower by a factor of ~10 times than those found for stars of similar 



- 2 - 



mass on convective tracks. The dramatic decrease in flux from convective to 
radiative phases of PMS evolution is likely related to major structural changes 
which influence the efficiency of magnetic field generation and thus the level of 
magnetic activity. 

As in the ONC, we find that stars with measured periods have, on average, 
higher X-ray luminosities. However, there is a wide range in and L^/ Ly^^i for 
both periodic and non-periodic stars of similar mass. Among stars with measured 
periods, the level of X-ray emission does not correlate with the rotation rate. 

For this data set, we find no statistically significant correlation between X-ray 
flux and (a) the presence or absence of circumstellar accretion disks, or (b) disk 
accretion rates as assessed from ultraviolet excesses. 

Subject headings: (stars:) (Galaxy:) open clusters and associations: individual 
(NGC 2264, Orion) 

1. Introduction 

Low-mass stars of all stages of evolution — from protostars to the main sequence — emit 
X-rays (see, e.g., Feigelson & Montmerle 1999). The relationships between the levels of X-ray 
emission and stellar properties for pre-main-sequence (PMS) stars are, however, very different 
from the well established correlations between rotation, Rossby number, and X-ray emission 
seen in young main-sequence stars (e.g., Stauffer et al. 1994, Micela et al. 1996, Pizzolato 
et al. 2003). Several factors may contribute to these differences. Low-mass PMS stars are 
fully convective during their early evolution, and differences in stellar structure (convective 
in the PMS phase vs. radiative core, convective envelope on the ZAMS) are likely to produce 
different levels of magnetic-field-related activity. Many PMS stars are still surrounded by 
accretion disks, which may influence X-ray emission in at least three ways. Additional X-ray 
emmission above coronal levels may be produced in accretion shocks or wind/jet shocks (e.g., 
Kastner et al. 2002, 2005; Guedel et al. 2005). Depending on the system orientation and 
the size of disk (and extent of flaring), disks or neutral winds may absorb X-ray photons, 
decreasing the detected flux, particularly at short wavelengths (e.g., Kastner et al. 2005; 
Walter & Kuhi 1981). The rotation speed of PMS stars could, in principle, also affect X-ray 
emission. If so, then the X-ray properties of stars surrounded by accretion disks may differ 
systematically from their counterparts that lack disks, owing to the effects of "disk locking" 
on rotation {e.g., Shu et al. 1987, 2000; Konigl 1991, 1989). Finally, there might also be 
other less direct interactions between accretion, the young stellar object, its outer convective 
zone, the magnetic field structures, and the X-ray emitting plasma, as discussed in Preibisch 
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et al. (2005). 

The recent COUP project (Getman et al. 2005) has investigated the detailed relation- 
ships among stellar parameters and X-ray flux among solar-like PMS stars in the ONC 
(Preibisch et al. 2005), which has an age of ~0.5 Myr. We consider here the relationships 
among rotation, X-rays, and disks for several hundred stars in the Orion Flanking Fields 
(~1 Myr) and NGC 2264 (~3 Myr). We focus our analysis on objects of masses ~0.3 to 2.5 

Mq. 

The stars in these regions are, on average, slightly older than those in the ONC and are 
therefore, in combination with the COUP data, well-suited for searching for changes during 
the early evolution of PMS stars. Importantly, for stars of ^^1-2 Mq, the current sample 
includes stars spanning the range of ages in which substantial structural changes take place 
as stars evolve from a fully convective phase to a phase in which they develop radiative 
cores and convective envelopes. Since the origin of the X-rays in the Sun is believed to be 
linked to dynamo-driven activity, we might expect that the mechanism of magnetic fleld 
and consequent coronal X-ray generation changes as stars make the transition from a fully 
convective state to one in which a radiative core is surrounded by a thin convective envelope. 
Consequently, we might also expect to see a change in the X-ray properties. 

The basic presentation of the data for the Orion FF and NGC 2264 appears in other 
papers (Ramirez et al. 2004a,b; Flaccomio et al. 2005a,b); we discuss here the imphcations 
of these data. After a brief summary of the observations (§2), we consider in some detail 
the completeness and selection effects of the sample of stars to be discussed here. Then we 
compare the relationships of log and log Lx/Lboi (§3) with mass and age for stars on 
convective tracks in all three clusters and find that the relationships are entirely consistent, 
thereby indicating that there is little or no evolution in X-ray properties during the first ~5 
Myr of evolution down convective tracks. Wc find clear evidence for a decrease in log Lx/Lboi 
when stars make the transition from convective to radiative tracks (§4). We also discuss the 
sample with rotation information (§5), and find that the subset of stars with measured 
periods has on average significantly higher values of log Lx/Lboi, as was found previously for 
the ONC (Preibish et al. 2005b). Also in accordance with previous studies, we find no clear 
trends of X-ray emission with rotation among the sample having spot-modulated periods. 
Finally, in §6, we examine the relation between disks, mass accretion rate (inferred from 
ultraviolet excess emission) and X-ray luminosity; we find no clear trends. 
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2. Observations, Data, and Sample Selection 
2.1. Observations 

In this paper, we discuss stars detected in four Chandra/ACIS fields: two in NGC 2264 
and two in the outer Orion Nebula Region, the so-called "Orion Flanking Fields" (FF; see 
RebuU et al. 2000, RebuU 2001 for more discussion of these fields). 

The Chandra data for the Orion FF were first presented in Ramirez et al. (2004b). The 
two fields here correspond to portions of fields 2 and 4 from RebuU et al. (2000); they are 
centered on 05h35ml9s, -4°48'15" and 05h35m6s, -5°40'48". Standard CIAO procedures 
were used to reduce the data. With about 48 ks exposure time per field, 417 sources were 
detected in total over both fields. The distance modulus we assumed for Orion is 8.36 (470 
pc). 

The Chandra data for the northern field in NGC 2264 were first presented in Ramirez 
et al. (2004a), in which 263 sources were detected in 48.1 ks. The ACIS-I array was centered 
at 6h40m48s, +9°51'. These data were also independently analyzed by Sung et al. (2004) in 
an effort to identify X-ray emitting WTTS. The ACIS data for the southern field in NGC 2264 
is discussed in Flaccomio et al. (2005, 2006). The ACIS-I array was centered at 6h40m58.1s, 
-|-9°34'00.40", and 420 sources were detected in 96 ks. Standard CIAO procedures were used 
to reduce the data. The distance modulus we assumed for NGC 2264 is 9.40 (760 pc). 

Following the approach discussed in Ramirez et al. (2004a,b) and Flaccomio et al. (2006), 
we determined upper limits for all stars with optical data in these fields. 

ACIS data have also been recently published in the ONC core region by Feigelson et al. 
(2003b) and as part of the COUP project (Getman et al. 2005, Preibisch et al. 2005b). Those 
much deeper (~ 20 x) observations of an extremely rich and slightly younger region (Ramirez 
et al. 2004b) provide a dataset of great statistical weight with which we can compare our 
results. 



2.2. Ancillary Data and Counterparts 

In order to interpret these X-ray data in the context of young stars, we assembled 
extensive catalogs of mult i- wavelength photometric and spectroscopic data for each cluster, 
as discussed in greater detail in Ramirez et al. (2004a, b). These catalogs are hereafter 
referred to as the ancillary data catalogs (ADC). For NGC 2264, the bulk of the data in the 
ADC comes from RebuU et al. (2002), as weU as Park et al. (2000) and the Two Micron AU 
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Sky Survey (2MASS). The ADC also contains preliminary proper motions from B. Jones; see 
Rebull et al. (2002). We have added to the ADC recently published data from Lamm et al. 
(2004), Sung et al. (2004), and Dahm & Simon (2005). For the Orion Flanking Fields, the 
ADC was assembled from ~30 published articles ranging from the U magnitudes in Rebull 
et al. (2000) to the JHK photometry from 2MASS; see Ramirez et al. (2004b). The ADC 
has also been updated with more recently published photometric data from Sicilia-Aguilar 
et al. (2005), as well as with 126 additional vsini measurements and upper limits currently 
under preparation for publication (Rebull et al. 2006 in prep). 

For each of the clusters, wc matched each X-ray detection from Ramirez ct al. (2004a,b) 
to optical and/or NIR counterparts from the ADC based on comparison of RA and Dec of 
the individual objects combined with manual inspection of the near-IR field around each 
X-ray source. 

Following the same approach as in, e.g., Rebull et al. (2000), using the stars' positions 
in the dereddened I /V — I color-magnitude diagram (CMD), we calculated Lboi following 
the precepts summarized by Hillenbrand (1997). Similarly, we assigned masses to the stars 
according to the two different sets of PMS stellar models: (1) those of D'Antona & MazziteUi 
(1994, 1998; DAM), with conversion from the theoretical to the observational plane (includ- 
ing calculation of Lboi) using transformations found in Hillenbrand (1997); and (2) those of 
Sicss ct al. (2000; SDF), the Z=0.02 model with no overshooting, using conversions from 
Kcnyon & Hartmann (1995). The DAM models arc what we have used in previous papers 
(e.g., Rebull ct al. 2000), but the SDF models were used by COUP (Prcibisch ct al. 2005b); 
masses from both models are discussed here to facilitate comparison with these and other 
papers. If stars fell outside the grid provided by the models, no masses or ages were assigned; 
such stars were excluded from this analysis. Additional discussion regarding uncertainties 
in masses and ages appears in Rebull (2001); we repeat here that the largest uncertainties 
are in the models themselves, and the ages and masses used here are perhaps best thought 
of as proxies for relative placements in the observational CMD. 

In Ramirez et al. (2004b), we have found that the X-ray selected sample in Orion FF 
(~1 Myr) is on average slightly older than that for the ONC (~0.5 Myr), but not as old as 
that for NGC 2264 (~3 Myr). Repeating this analysis now including the second NGC 2264 
field, it is still the case that the median age of the FF sample is older than the ONC, and 
NGC 2264 is older still. The range of ages found within each of these clusters is on the order 
of ~3 Myr and comparable in all three clusters, with the largest fraction of older stars being 
found in NGC 2264; see Figure 1. The ages listed here (0.5, 1, and 3 Myr) are obtained from 
comparing distributions of ages derived using DAM models for the "best possible sample" 
(see next subsection below) for these clusters. If we use SDF models for the same subset of 
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Fig. 1. — Color- magnitude diagram (CMD) for the final best sample to be considered (see 
text). Dots indicate detections, and triangles indicate Lx upper limits. Isochrone lines 
overplotted are, from top to bottom, isochrones for 1 Myr and 3 Myr, and the ZAMS; the 
left panel is DAM98 and the right panel is SDF. Lines of constant mass are, left to right, 
2, 1, 0.5, and 0.15 M©. The apparent quantization arises from correcting the V-I colors to 
those expected from the spectral type. The range of ages found in these clusters is ~3 Myr, 
but the largest fraction of older stars is found in NGC 2264. The red symbols denote those 
stars that are no longer on convective tracks; see discussion in §4 below. 

stars, we find ages that are on average older by ~1 Myr, but the relative age placement of 
the clusters is similar at 1.5, 2, and 4 Myr. The relative placement of the clusters is also 
consistent with the relative numbers of radiative and convective stars (see below); the fewest 
radiative stars are found in the ONC and the most are found in NGC 2264. 



2.3. Catalog Statistics and Sample to be Considered Here 

In order to interpret our X-ray data, we have to include data from the ancillary data 
catalog (ADC). However, we cannot simply take all of the ADC data merged with the X-ray 
data and proceed. Not all stars have a complete set of data in the ADC, and there is a wide 
range of ADC survey depths among our sample clusters. In this section, we describe what 
kinds of data we require for a given star to appear in subsequent analysis in this paper. In 
the end, we settle on a sample which we will hereafter refer to as the "best possible sample" 
(as defined below and summarized in Table 1). Within this best possible sample, there are 
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stars that are on radiative tracks according to either the DAM or SDF models, and stars 
that are on convective tracks according to those models. These samples are referred to as the 
best possible radiative and best possible convective samples. In plots that follow, where both 
the radiative and convective stars appear, the radiative stars are colored red; otherwise, the 
plots contain just the convective stars. The following discussion is summarized in Table 1, 
where the numbers of detections and/or upper limits are tabulated for each of a variety of 
criteria. 

While there arc 416 X-ray sources detected in the Orion FF and 683 sources detected 
in NGC 2264 (hnc A in Table 1), only 84% (348) of the FF sources and 75% (509) of the 
NGC 2264 sources have optical and/or NIR counterparts in the ADC (line B in Table 1). 
For determining, among other things, masses and ages, we require a V and / magnitude so 
that we can place all the stars in the same CMD and at least derive self-consistent values. 
There are 278 (67%) of the FF X-ray detections and 465 (68%) of the NGC 2264 X-ray 
detections meeting those criteria (line D in Table 1). Further, of those, 94% (261 in the FF 
and 438 in NGC 2264) are likely cluster members based on position in the sky and position 
in the CMD (line E in Table 1) — these objects appear in a clear "locus" above the ZAMS 
(see RebuU et al. 2000, 2002 for full CMDs and additional discussion). 

Of the entire set of X-ray detections, only 115 (27%) Orion FF sources and 294 (52%) 
NGC 2264 sources can be identified with counterparts that have known spectral types in the 
ADC (line C in Table 1). Because of the original motivations behind the studies amassing 
most of these spectral types, the set of stars with types is strongly biased toward the stars 
with known rotational periods. As we will see below, stars with measured periods are brighter 
in X-rays, which therefore affects the set of stars with both X-ray detections and spectral 
types. 

As a result of this bias, we chose not to limit our analysis here to solely the stars 
with spectral types. In order to minimize the chances of field star contamination, we have, 
however, restricted our discussion to objects (261 detections in the FF and 438 detections in 
NGC 2264) that are likely members based on position in the sky and in the CMD (line E in 
Table 1). 

Upper limits can be important for determining trends in complex data sets such as this. 
We obtained upper hmits in Lx for stars with optical detections in the ADC. For the FF, 
this adds 84 upper hmits to the 261 detections to be considered (for just the stars that are 
hkely to be members - the ones in the CMD locus - an increase of ~30%). The optical 
surveys of NGC 2264 incorporated into the ADC are very deep and are not well-matched 
to the relatively shallow Chandra survey. If the entire optical sample from the ADC were 
included (line F in Table 1), the number of X-ray upper limits would dwarf the detections 
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(and moreover would likely be dominated by non-members). We have, therefore, limited our 
analysis to likely member stars (line G in Table 1; those both on the molecular cloud and 
in the CMD locus) with / < 16, which is well-matched both to the Lx detections and to 
the optical surveys incorporated into the ADC for both clusters. In order to make a fair 
comparison with the set of detections, we further limit even the detections to those stars 
with / < 16. This adds 168 upper limits to the now-reduced 349 detections to be considered 
in NGC 2264 (line H in Table 1). These restrictions should also reduce the contamination 
of our sample by non-members. The percentage of non-members is probably higher among 
stars with no detections (because members are more likely to be bright in X-rays) and at 
the faint end (because the mass function means there are simply more low-mass stars in the 
field). This cut in I of course has implications for the mass and age range to be considered 
here; see Figure 1. This cut eliminates stars with M <0.3 Mq and t >5 Myr, using DAM 
models. 

To summarize the discussion to this point (line H in Table 1), then, there are 261 X-ray 
detections with optical counterparts from the ADC that are likely stellar members of the 
Orion FF, with 84 additional upper limits. There are 349 X-ray detections with optical 
counterparts from the ADC with / < 16 that are likely stellar members of NGC 2264, with 
168 additional upper limits. 

Because we will be comparing the X-ray properties of stars on radiative tracks with 
the stars on convective tracks, we have further limited the discussion to just stars that fall 
within the range of L and Teg spanned by published stellar models (the "model grid"), and 
those with M <2.5 M©. This leaves 250 detections and 80 upper limits in the FF, and 317 
detections with 139 upper limits in NGC 2264 (lines I and J in Table 1) 

The sample considered here uses the best available information for each star, e.g., if 
there is a spectral type in the ADC, the star is dereddened by values specific to each star, 
and if there is no type available, the star is dereddened by the most likely reddening in 
the direction of the cluster (see RebuU et al. 2000, 2002 for more discussion on reddening 
towards these clusters). A CMD for this sample appears in Figure 1. This constitutes the 
"best possible sample" for both clusters. Note that there are no stars included in this sample 
that fall outside the model grid. Note also that despite using two different models, DAM 
and SDF, the number of stars in the best possible sample using DAM models is very similar 
to the best possible sample using SDF models. 

For each of the evolutionary models (DAM & SDF) depicted in Figure 1, the red symbols 
indicate the stars between ~1 and ~2.5 Mq that are on radiative tracks (lines K and L in 
Table 1). There are very few of these radiative stars available in the FF, with only 1 X-ray 
detection (and 1 upper limit) using DAM models, and 18 X-ray detections (and 1 upper 
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limit) using SDF models. But, in NGC 2264, using the DAM models, there are 17 X-ray 
detections and 16 upper limits; using the SDF models, there are 62 X-ray detections and 24 
upper limits . The red, radiative points here are the "best possible radiative sample. " 

Essentially all of the discussion that follows is restricted to stars on convective tracks 
with masses <2 M©, e.g., the black points in Figure 1. The black, convective points here 
are the "best possible convective sample. " For the best possible convective sample, using 
DAM models, this leaves 249 detections and 79 upper limits in Orion so 24% of the best 
possible convective sample is upper limits; in NGC 2264, there are 300 detections and 123 
upper limits, so 29% of the best possible convective sample in this cluster is upper limits. 
Similar numbers are obtained for SDF models; see Table 1. In most dicussions that follow, 
we consider only the convective sample. Where the radiative stars are included in subsequent 
discussion and plots, we note this explicitly. In plots that follow, where both the radiative 
and convective stars appear, the radiative stars are colored red. 

Of this final, best possible sample of convective stars, we also want to know how many 
stars have measured periods and vsini values (fines O and P in Table 1). There are 110 
stars in the Orion FF with measured periods that are also detected in X-rays, and 14 stars 
with measured periods but upper limits in X-rays; there are 64 stars with measured vsini 
and X-ray detections, and 4 with upper limits in X-rays. In NGC 2264, there are 181 stars 
with measured periods and X-ray fluxes, and 38 more with upper limits in L^. There are 
only 30 stars with measured v sin i and X-ray detections, and 3 more stars with upper limits 
in Lx. 



Table 1. Summary of selection criteria used to assemble "best possible sample. 
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3. Lx vs. Lboi 

All recent investigations, including the COUP study (e.g., Preibisch et al. 2005b), have 
found a clear positive correlation between log and log Lboi but with a very large intrinsic 
(as opposed to instrumental) scatter (Feigelson et al. 2003). 

Figure 2 shows the results for the ONC from the COUP sample, Orion FF, and NGC 
2264. The best possible radiative sample for NGC 2264 and the Orion FF is shown in red, 
and the best possible convective sample is shown in black (essentially all of the COUP stars 
are convective). Also shown is the best fit straight line obtained by Preibisch et al. (2005b) 
for the ONC COUP data, which has a sensitivity hmit of about log Lx=27.0. This should 
be compared with a sensitivity limit of about log Lx ~ 28.5 for our observations. The line 
derived for the ONC is a reasonable fit to the data for stars on convective tracks for both 
of our regions. We note, however, that many of our points near log Lboi~33.5 lie above the 
mean line for the ONC, and there is a cloud of points with low values of log that lie below 
the line near log Lboi~32.8, suggesting that perhaps a linear fit may not be appropriate. 
However, both of these features are also seen in the ONC plots and may thus refiect some 
inadequacy in the ability of a straight line to fit the data for both data sets. Figure 3 depicts 
the same data as in Figure 2, but in box plot form. These box plots have been used in 
other papers (e.g., Flaccomio et al. 2003c) as a mechanism for interpreting scatter plots. We 
emphasize that the boxes in Fig. 3 have been corrected using the Kaplan-Meier estimator 
for censored data to take into account the upper limits present in the data. The central 
line in each box denotes the (corrected) median; the ends of each box are the first and third 
quartile of the (corrected) distribution; the lines extend to the most extreme values that are 
not more than 1.5 times the interquartile range; and the open circles are those points outside 
1.5 times the interquartile range. If downward-pointing triangles appear in the boxes, then 
the true (corrected) lower limit of the box is located at some unknown location below the 
line indicated. The top three panels in Fig. 3 contain only convective stars; the bottom panel 
contains only the radiative stars in NGC 2264 (there are not enough radiative stars in Orion 
to merit a separate plot). Examination of Fig. 3 reveals that the X-ray properties of the 
radiative stars are significantly different than those of the convective stars. 

Preibisch & Feigelson (2005) used the Orion data to search for changes in X-ray emission 
with age. They assumed that the stars in Orion ranged in age from 0.1 to 10 Myr (Palla 
& Stabler 1999). Because contraction down the Hayashi tracks slows with increasing age, 
at ages of ~3 Myr. typical uncertainties in derived luminosities and effective temperatures 
lead to age uncertainties of ~2-3 Myr (lo" errors, Hartmann 2001). Based on the ONC 
data alone, Preibisch & Feigelson concluded that the mean decayed slowly with age (Lx 
(X r^^/^) during the first 10 Myr for stars with masses in the range 0.5-1.2 M©. However, in 
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Fig. 2. — Plot of log La; vs. logLboi for the stars (see §2.3) in the Orion ONC from the COUP 
sample (top plot), Orion FF (middle plot) and NGC 2264 (bottom plot). Only detections 
from COUP are shown, whereas upper limits are shown for the other two panels. The upper 
limits occur at all values of because the PSF varies significantly across the field of view, 
so a star appearing at the edge of the field will have a much higher upper limit than a star 
near the center of the field. The boundaries for this plot were optimized for the FF and NGC 
2264 data; 24 detections and 14 limits from the COUP data appear outside of this region, 
all but 4 detections (and 1 limit) of which appear below log Lx=28. The solid line is the 
best fit value to the ONC COUP data analyzed by Preibisch et al. 2005b. The dotted line 
is the canonical "saturation" value of logL^^/Lboi ~ —2.9. In the lower two panels, we see 
a suggestion of a bifurcation in the distribution beginning at about half the luminosity of 
the Sun and extending to higher luminosities. In NGC 2264, 50% of the stars on the lower 
branch are on radiative tracks (indicated in red) and therefore may have X-ray properties 
that are different from stars on convective tracks. A discussion of the comparison with the 
ONC is provided in the text. 
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Fig. 3. — Plot of log VS. logLboi for the stars from Figure 2, but represented in box plot 
form. The first three panels are convective stars; the last panel is the radiative stars in NGC 
2264 (there are not enough radiative stars in the other clusters to merit a separate plot). It 
is readily apparent here that the X-ray properties of the radiative stars are different than 
those of the convective stars. 
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order to match the X-ray properties of ZAMS stars in young clusters, the decay of must 
be accelerated at ages t >10 Myr. For stars with masses in the range 0.1-0.4 Mq they found 
only a modest decrease in log during the first 100 Myr and a possible slight increase in 
log Lx/Lboi over that same time period. Flaccomio et al. (2003c) also investigate trends of 
X-ray emission with time, comparing median levels of X-ray emission in several star-forming 
regions. They found constant log and an increase in log Lx/Lboi with time through ~3-4 
Myr and a subsequent leveling off at saturation level (log Lx/Lboi ~ —3). 

Our own data are for clusters which contain stars that arc on average older than the 
stars in the ONC. We see two differences with respect to the ONC data. First, the upper 
envelope in log Lx for ONC data Ues about 1.5 dex or a factor of 30 above the best fit straight 
line. For the Orion FF and NGC 2264, the upper envelope is only about a factor of 10 (1 
dex) above the best fit straight line. As Fig. 2 also shows, very few of the Orion FF and 
NGC 2264 stars he above the nominal "saturation" hmit of log Lx/Lboi = —2.9, while a large 
number of ONC stars lie above this limit (Preibisch et al. 2005b). Quantitatively, 14% of the 
detected ONC stars have log Lx/Lboi > —2.9, and 4.6% have log Lx/Lboi > —2.5. For the 
Orion FF, the fractions are 5% and 0.4%, respectively, and for NGC 2264, the fractions are 
7% and 0.9% respectively. In other words, there is a significantly larger number of stars with 
high Lx and high Lx/Lboi in the younger ONC. There are three possible explanations for 
this effect. First, it could be entirely an age effect, in that the most extreme levels of activity 
may decay during early PMS evolution. Second, it could be the intrinsic variability of the 
stars themselves in Lx. Third it could be the intrinsic variabihty of the stars themselves 
in the optical photometry that goes into the Lboi estimate. In either of these latter cases, 
the variability could be more extreme in the younger stars, allowing for more outliers. With 
regards to the variability in Lx, Favata et al. (2005), in studying the brightest ~1% of the 
COUP flares, flnd some objects that flare 1-2 orders of magnitude, and that the characteristic 
times for these flares range from 10-400 ks. In the spcciflc case of the COUP sample, however, 
as Preibisch & Feigelson (2005) point out, the influence of individual flares on the total Lx 
reported for any one star is reduced by the long time baseline of their observations. Thus, 
outliers solely due to flares in the COUP data set are less likely than in the other (much 
shorter) observations. Variability in the optical data that go into calculation of Lboi may 
account for some fraction of the outliers; optical variability can exceed a magnitude. RebuU 
et al. (2002b) discussed the sources of these errors in detail for stars with spectral types over 
limited mass ranges and in regions of less reddening. For the COUP sample in particular, 
even though reddening is a concern, spectral types are known for every object, and so it 
is unlikely that this can account for all (or even most) of the outliers. More data on time 
variability will be required to determine why there are more outliers with high Lx in the 
ONC. 
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The lower envelope of the X-ray detected stars in all three regions for both Lx and 
Lx/Lboi, however, is comparable for stars with Lboi similar to that of the Sun, and so the 
decay in the mean and lower bounds of the X-ray luminosities with time during the first 5 Myr 
of PMS evolution indeed appears to be small, as Preibisch & Feigelson (2005) reported for 
the ONC alone. Unfortunately, because of the shallow depth of our current survey, we cannot 
compare the lower bounds of the detections of stars with M < 1 M© in these three regions. 
Since more than half the stars in our sample were detected, we can compare the median 
log Lx and log Lx/Lboi for the best possible convective sample in all three clusters. The 
median log for the ONC (further restricted to M >0.16 Moto account for the difference 
in survey depths), the Orion FF, and NGC 2264 are 29.55, 29.71, and 29.90, respectively; 
the median values of log Lx/Lboi (in the same order) are —3.71, —3.50, and —3.44. (There 
are some upper limits in this sample, so as a simple worst-case test, we recomputed the 
medians, adopting the upper limits for all stars as detections, and it does not change these 
median values very much.) These results are consistent with that obtained by COUP for the 
ONC. While it is true that the most luminous objects have disappeared by the age of NGC 
2264, there is not much change in the median values. We certainly do not see a significant 
decrease during the first 5 Myr, in agreement with the Preibisch & Feigelson (2005) work on 
the ONC alone. 



4. X-ray properties and stellar structure: the transition from convective to 

radiative tracks 

In this section, we compare the X-ray properties of low mass stars still located on fully 
convective tracks with those of higher mass stars, some of which appear to lie on convective 
tracks, while others of slightly greater age have developed radiative cores and are proceeding 
toward the main sequence along radiative tracks. In Figure 1, the stars on radiative tracks 
are indicated. 

In Figures 4 and 5, we plot log Lx and log Lx/Lboi vs. mass, retaining the red (radiative) 
and black (convective) convention. Among the detections, the stars on radiative tracks have 
lower values than their similar mass counterparts located on convective tracks; moreover, 
there is a higher fraction of upper limits among the radiative stars than among the convective 
stars. Examining solely stars with detections in NGC 2264, the median log Lx/Lboi for 
all convective stars between 1 and 2.5 is —3.51; the median log Lx/Lboi for radiative 
stars (over the same mass range) is —4.01. However, many of the values for radiative stars 
are upper limits; adopting the upper limits for all stars as detections at the level of their 
upper limits, the median log Lx/Lboi for the convective stars changes subtly to —3.54, while 
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Fig. 4. — Plot of log La; (top plots) and logL^/ L^oi (bottom plots) vs. log M/M© from DAM 
for the Orion FF (left plots) and NGC 2264 (right plots), for the best samples available for 
each cluster (see §2.3). The solid line is the relationship between and SDF masses found 
for the ONC by Preibisch et al. 2005b. Note again the apparent bifurcation at about 0.5 
Mq. Red symbols are those stars still on radiative tracks; see Fig. 1. Dotted lines indicate 
0.1, 0.16, 1, and 2.5 Mg, from left to right. 
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Fig. 5. — Plot of logLrc (top plots) and logL^^/Lboi (bottom plots) vs. log M/M0 from SDF 
for the Orion FF (left plots) and NGC 2264 (right plots), for the best samples available for 
each cluster (see §2.3). The solid line is the relationship between and SDF masses found 
for the ONC by Preibisch et al. (2005). Note again the apparent bifurcation at about 0.5 
Mq. Red symbols are those stars still on radiative tracks; see Fig. 1. Dotted lines indicate 
0.1, 0.16, 1, and 2.5 Mq, from left to right. 
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the median for the radiative stars changes to —4.48. By assigning the upper hmits to be 
"detections," we of course obtain a conservative hmit on the median log Lx/Lboi; even so, 
the medians for the convective and radiative samples differ by a factor of 10. Evidently the 
transition from a fully convective state to a radiative-core convective envelope configuration 
results in a dramatic decrease in coronal activity as diagnosed by X-rays; Gagne et al. (1995) 
also report differences in the X-ray properties between the convective and radiative stars in 
Orion based on ROSAT observations (see also Strom 1994 and Flaccomio et al. 2003b). 

We have examined the COUP data to search for a similar trend among the ONC sample. 
Unfortunately, because the region is younger, too few stars have evolved to their radiative 
phases to enable a convincing test. 

Preibisch et al. (2005b) reported a weak but significant correlation between log Lx/Lboi 
and mass. Spearman's correlation coefficient does not find any significant correlations in our 
data for the Orion FF and NGC 2264. However, we note that our survey covers a much 
smaller mass range than COUP. Restricting the comparison with COUP data to the same 
mass range reveals no clear trends in either sample. Because our data are dominated at low 
masses by upper limits, we are also unable to search for mass-dependent differences in log 
Lx/L\,oi as reported by Flaccomio et al. (2003b, 2003c), who find similar log Lx vs. mass 
relationships to that reported by Preibisch et al. (2005b). 

The conclusion that X-ray luminosity is related to the interior structure of PMS stars is 
of course vahd if and only if stars assigned to radiative tracks are placed correctly in the HR 
diagram, are predominantly members of the Orion FF and the ONC, and if the X-ray flux 
is not dominated by a faint companion that is not visible in the optical. We now examine 
each of these requirements for NGC 2264, which has the largest number of massive stars in 
apparent post-convective evolution. 

How accurately have we located these stars in the HR diagram? These stars are among 
the brightest stars in the cluster, and for 94% of the radiative stars (both detections and 
upper limits), we have spectral types; they range in type from AO to K4. Because we have 
types, we are able to deredden each star individually with a correction appropriate to that 
specific star. Therefore, the placement of these stars in the dereddened CMD would appear 
to be as secure as possible for young (often variable) stars. These stars appear in the correct 
vertical position to be taken as cluster members; they do not appear significantly above or 
below the ZAMS for a cluster at the distance of NGC 2264. We therefore believe that most 
of these stars are members based on position in the HRD. 

Because these stars are among the brighter stars in the cluster, we have preliminary 
proper motions (see RebuU et al. 2002) for 54% of the stars on radiative tracks; 33% of 
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the red points (detections and upper limits) in Figure 1 (DAM models) are likely members, 
22% are likely non-members, and the rest have no proper motion membership information 
available. Of the available sample, at least 60% of the stars are thus members based on 
proper motions. 

Another indicator of likely membership is the presence of a circumstellar disk. Seventy 
percent of the best possible radiative sample has positive circumstellar and/or accretion disk 
indicators, suggesting at least 70% of the sample are indeed members. 

Are unresolved binaries important? They may play a role, particularly in the case of 
X-rays where low-mass companions can outshine their primaries. However, what we have 
shown is that the fluxes of the stars on radiative tracks are systematically lower than those of 
stars with similar masses on convective tracks. If optically fainter companions do contribute 
some of the flux, then the discrepancy must be even larger. 

In summary, then, we believe that most (>70%) of the best possible radiative sample 
are indeed cluster members. We therefore conclude that there is a significant decrease in 
X-ray luminosity as stars of similar mass evolve from convective tracks, where log L^/ L-^oi 
values lie close to the "saturation limit," to radiative tracks, where L^/L\^o\ values arc ~10 
times smaller. In support of this result we note two additional observations: 1) the apparent 
'bifurcation' in Lx/Lboi disappears for masses lower than about 0.5 M© which is approxi- 
mately the mass below which stars remain fully convective from the stellar birthline to the 
ZAMS; and 2) the levels of Ly_ and Lx/Lboi for stars on radiative tracks are roughly in agree- 
ment with the values found in this mass range in the NEXXUS survey of nearby field stars 
(Schmitt & Liefke 2004). 



Table 2. Summary of recent studies of large surveys of star- forming regions considering effects of 
disks/ accretion/degree of embeddedness or period on Lx 



work 


instrument 


cluster 


range of L^/Lj^^i'^ 


disk indicator 


conclusion on disks^ 


conclusion on rotation 


Bouvier (1990) 


Einstein 


Taurus- ^Vuriga 


logLx=27-32 


Ha >10A 


no difF. 


anti"Correl. 


Feigelson et al. (1993) 


ROSAT/PSPC 


Cham I 


—2 to —5 


EW(Ha)>10A 


Lx{C)<Lx(yV) but no diff. within sel. efF. 


no correl. 


Gagne &c Caillault (1994) 


Einstein 


ONC 


—2 to —7 


H — K, Ha 


no difference 


no correl. 


Damiani Sz Miccla (1995) 


Einstein 


Taurus- Auriga 


(—2 to —4.5) 
—3 to —5 


IRAS-25 


WTTS fainter 


weak anti-correl. with scatter 


Gagne et al. (1995) 


ROSAT/HRI 


ONC 


—2 to —7 


A(H — K) > 0.1 


Lx, Lx/Lbol(C)<La:, Lx/Lbol(y^) 


no correl. 








(—2 to —5) 








Casanova et al. (1995) 


ROSAT/PSPC 


p Oph 


—2 to —5 


Class 0-3 


(®™^bedded) "^Z/jj {less embedded), no diff. 


N/A 


Neuhauser et al. (1995) 


RASS 


Taurus- Auriga 


-4 to -7 


EW(Ha)> 5-15 A 




correl. 












I/x(C) harder than W 




Lawson et al. (1996) 


ROSAT/PSPC 


Cham I 


—2 to —5 


Ha 


no difference 


N/A 


Preibisch (1997) 


ROSAT/PSPC 


several 


logi:,x=27-32 


N/A 


N/A 


anti-corrcl. 


Alcala et al. (2000) 


RASS 


Orion 


logLx~31 


N/A 


N/A 


no correl. for K stars 


Wichmann et al. (2000) 


RASS 


Taurus- Auriga 


—2.5 to —5 


N/A 


N/A 




Flaccomio et al. (2000) 


ROSAT/HRI 


NGC 2264 


—2 to —5 


Ha 


no difference, but C more variable 


N/A 








(—2 to —4) 








Grosso et al. (2000) 


ROSAT/HRI 


p Oph 


—2 to —4 


ISO data 


no difference 


N/A 


Stclzer et al. (2000) 


ROSAT/PSPC 


Taurus- Auriga 




HBC listing W/C 


Lx(G)<Lx(W) but stronger flares, 


anti-correl. 












Sz C more variable 




Imanishi et al. (2001) 


Chandra/ACIS 


p Oph 




Class 0-3 


Lj^/L 1,^1 (Class I)<L^/Lboz(Class II & III) 


N/A 


Stelzer & Neuhauser (2001) 


ROSAT/PSPC 


Taurus- Auriga 


-3 to -6 


EW(Ha)> lOA 


i.x(C)<Lx(W) 


correl. 


Preibisch & Zinnecker (2002) 


Chandra/ACIS 


IC 348 


-2 to -6 


EW(Ha)> 10 A, 


Z,x(C)<Lx(W), but likely to be sel. eff. 


N/A 








(-2 to -4.5) 


A(K - L) > 0.17 


(goes away for A.(K — L) 




Feigelson et al. (2002, 2003) 


Chandra/ACIS 


Orion (ONC) 


-2.5 to -6 


A(J - K) > 0.3 


no difference 


weak correl. 


Getman et al. (2002) 


Chandra/ACIS 


NGC 1333 




JHK 


no difference 


N/A 


Tsujimoto ot al. (2002) 


C]handra/ACIS 


OMC-2, OMC-3 




Class 0-3 


fcT of C higher than W, but poss. scl. off. 


N/A 


Flaccomio et al. (2003a) 


ROSAT/HRI 


NGC 2264, 


-2 to -7.8 


A(/ - K) > 0.8 


L:,,L^/Li,^l(C)<L^,L^/Li,,l{W) 


N/A 






Cham I 


(-2 to -5) 








Flaccomio et al. (2003b) 


Chandra/HRC 


ONC 


-2 to -7.8 


EW(Ca II) > 1 A 


La:,L^/Lt,^l(C)<L^,L^/Li,^l{W) 


no correl.; stars with P high 


Stassun et al. (2004) 


Chandra/ACIS 


ONC 


(-2 to -5) 
-2 to -6 


EW(Ca II) > 1 A 


L^,L,,/Li,^l{C)<L^,L,,/Li,^l(W) 


log Z>x/^bol correl. with P, but not 














vsini; stars with P high Lx 


Feigelson & Lawson (2004) 


Chandra/ACIS 


Cham I 


=28-31 


various 


Lj:(C)<L:,{W) (slightly) 


N/A 


Preibisch (2005) 


XMM/EPIC 


Serpens 


3 


literature 


Class I more variable than Class II + III 


N/A 


Ozawa et al. (2005) 


XMM/EPIC 


p Oph 


-2 to -5 


literature 


Class I higher T, absorp. than II+III 


N/A 


Preibisch et al. (2005b) 


Chandra/ACIS 


ONC (COUP) 


-1 to -6 


A(/ - K), 


accretors less active 


weak correl. 










A{K - L), EW(Ca II) 






^Values in parentheses are for 


the low-mass stars 


in the sample. 











'C" = CTTS; "W" = WTTS 
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Fig. 6. — The log Lx/Lboi distribution for the best samples available for each cluster (see §2.3, 
including only convective stars) for stars with and without measured periods. Stars with 
periods are X-ray bright, and stars that are X-ray bright are more likely to have measured 
periods. The line at an angle is a fit to the ratio. The slopes are 0.34±0.05 (Orion FF) and 
0.26±0.05 (NGC 2264). 



5. X-ray Emission and Rotation 

Flaccomio et al. (2003b) first found evidence that stars in the ONC for which rotation 
periods are known to have higher X-ray luminosities than stars which show no clear evidence 
of spot-modulated rotation; Stassun et al. (2004) confirm this result. 

The same trend is seen among fully convective stars (M < 0.5 Mq) in both the Orion FF 
and NGC 2264. Figure 6 shows that the stars with derived periods in these two associations 
are clearly brighter in log Lx/Lboi, and, conversely, stars that are bright in log Lx/Lboi are 
more likely to have a measurable period. For the best possible convective sample (plotted in 
Figure 6), the best-fit slopes to the lower panels in the Figure are (for the FF) 0.34±0.05, and 
(for NGC 2264) 0.26±0.05. These values do not change significantly when considering other 
possible subsamples of the data. Median values of log L^/L-^^i for the best possible convective 
sample for the ONC, FF, and NGC 2264, for the sample with periods, are (respectively) 
—3.57, —3.35, and —3.37. For the best possible convective sample without periods (in the 
same order), the medians are —3.81, —3.65, and —3.53. It is clear that the samples with 
measured P are brighter in X-rays (and again there are no trends with age). 
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log (velocity or vsini/[km/s]) 

Fig. 7. — Cumulative distributions of measured vsmi (multiplied by A/tt to account for 

projection effects; dotted line) and calculated velocity (solid line) for the FF (top panel) 
and NGC 2264 (bottom panel). Aside from where the spectral resolution cuts off the vsini 
measurements, the distributions are indistinguishable. 
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In searching for correlations between X-ray emission and rotation rates of PMS stars, 
another concern is whether the sample of stars with known periods represents the full range 
of rotation rates for PMS stars. It might be, for example, that periodic stars are biased 
toward PMS stars dominated by either long or short rotation periods, or by stars viewed 
more nearly equator-on than pole-on. The study of periods and v sin i in the ONC by Rhode 
et al. (2001) concluded that there was no difference in the distribution of rotation rates for 
samples of stars with periods and samples of stars with v sin i. A similar test is performed for 
the current data set in Figure 7. For the purpose of this graph, rotational velocities have been 
derived from P and an estimate of the radius (see, e.g., RebuU et al. 2002 for more on this 
calculation); the vsini values have been converted to velocity by multiplying the observed 
rates by 4/7r to allow for inclination effects (Chandrasckhar & Muench 1950; Gaige 1993). 
The Orion FF appears to exhibit a relative paucity of objects with low vsmi. However, this 
reflects the resolution limit of the data rather than a real physical effect. Taking this into 
account, we conclude from these data that there is no apparent bias in the period-selected 
sample. However, our data set is too sparse to distinguish whether the vsini distribution 
for the low Lx/I/boi stars that lack periods differ from the ensemble distribution. 

Why are we more likely to be able to determine periods from photometric monitoring 
programs for stars with large values of log Lx/Lboi? In order to obtain an optical photomet- 
ric period, we require a clear view of the stellar photosphere where there are large, stable, 
relatively isolated spots or spot groups rotating into and out of view causing periodic mod- 
ulations. Stassun et al. (2004) have suggested that perhaps more active stars have larger 
and/or more organized magnetic spot coverage and that periodic variations are more easily 
measured in such circumstances. If so, then stars with more organized spots also have higher 
levels of coronal activity leading to higher X-ray luminosity. 

Since the earliest days of X-ray astronomy, investigators have sought a relation between 
rotation rate (periods or vsvai) and X-ray luminosity for PMS stars. Because the X-ray 
flux is thought to originate in coronal activity which is driven by convection and rotation, 
and because faster-rotating main sequence stars are brighter in X-rays (e.g., Stauffer et al. 
1994, Micela et al. 1996, Pizzolato et al. 2003), it was perhaps reasonable to expect to find a 
similar relation in younger stars, but the evidence remains mixed. Table 2 contains a listing 
of the conclusions on this topic from the hterature to date. Most recently, the COUP project 
(Preibisch et al. 2005b) found a weak but statistically significant correlation between P and 
log Lx/Lboi- 

Figure 8 plots log and log Lx/Lboi against log period for our sample, and Figure 9 
shows the same data but in box plot form. There is no significant trend in these plots, nor in 
similar plots of X-ray emission against linear P, vsvai, or 1/P. Upper limits are found at all 
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Fig. 8. — Log Lx and log Lx/i^boi vs. log period for the best samples available for each cluster 
(see §2.3). There is no correlation here, but keep in mind that stars with periods are X-ray 
bright (see previous figure). The dashed line is the relationship between Lx/Lboi found by 
Preibisch et al. (2005) for the ONC, and the solid line is the approximate relationship found 
for main sequence stars. We do not see a similar correlation (either the ONC or MS one) for 
either the Orion FF or NGC 2264. 



periods and rotation rates. The current data do not support as strong a correlation between 
period and fractional X-ray luminosity as the one reported for the ONC by Preibisch et al. 
2005b. We also do not find the same correlations as can be found in main sequence stars. 



6. Disk and Accretion Effects 



Two processes associated with disk accretion can, in principle, affect the observed X- 
ray properties of solar-type PMS stars: (1) local heating of the stellar photosphere at the 
'footprints' of magnetospheric funnel fiows. where highly supersonic {v > 200 km sec^) 
gas lands on the stellar surface {e.g. Kastner et al. 2002; Schmitt et al. 2005), producing 
shocks, local heating and, as a result soft X-rays; and (2) absorption of keV coronal X-ray 
emission by gas associated with accretion-driven outflows (Walter & Kuhi, 1981). The flrst 
of these processes would qualitatively produce an increase in X-ray emission above coronal 
activity levels as disk accretion rates increase. Quantitative estimates of the magnitude of 
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log(Period/d) log(Period/d) 



Fig. 9. — The data from Fig. 8 but in box plot form. We do not see a correlation similar to 
that from the ONC (dashed line) or from the main sequence (sohd line) for either cluster. 
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Fig. 10. — Log Lx and log Lx/Lboi vs. I — K excess disk candidates for the best samples 
available for each cluster (see §2.3). Disk candidates are redder than 0.3 mags (to the right 
of the dotted line; see text). 
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Fig. 11. — Log Lx and log Lx/Lboi vs. U — V excess disk candidates for the best samples 
available for each cluster (see §2.3). Disk candidates are bluer than —0.5 mags (to the left 
of the dotted line; see text). 



the expected increase are not available, owing primarily to the complexity of the detailed 
radiative transfer studies needed to estimate the X-ray flux emerging from accretion spots. 
The second process would produce a decrease of X-ray emission with increasing accretion 
rate, since the column density, and hence the X-ray optical depth of outflowing material, is 
expected to increase with increasing accretion rate. For isotropic outflows, significant 5 keV 
absorption is expected to occur for mass outflow rates of 10~^ M© yr~^ or mass accretion 
rates of ~ 10"^ Mq yr'^ 

There has been a great deal of work attempting to uncover relationships between cir- 
cumstellar disks, and/or disk accretion, and/or X-ray emission (either log Lx or log Lx/Lboi)- 
Table 2 collects the conclusions from large surveys of X-ray emission for accreting and non- 
accreting PMS stars. Kastner et al. (2002); and Schmitt et al. 2005 find evidence of enhanced 
X-ray emission associated with transient increases in disk accretion rate, thus suggesting that 
at least some X-ray emission may be associated with accretion as opposed to coronal activity. 
In ensemble, there is no evidence of excess X-ray emission for accreting PMS stars. Rather, 
the trend, if any, is for an anti-correlation between the presence of an accretion disk and 
total X-ray luminosity. 

Walter & Kuhi (1981) were the first to report an inverse correlation between soft X-ray 
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flux and Ha emission for classical T Tauri stars. Working under the assumption that the 
observed equivalent width of Ha serves as a surrogate for mass outflow rate, they argued that 
the apparent decrease of X-ray emission with increasing Ha equivalent width reflected the 
effects of X-ray absorption by outflowing gas. More recently, Flaccomio et al. (2003b) find a 
similar result: a larger fraction of stars with active accretion signatures, there measured by 
Ca II emission, have fainter values than non-accretors when stars of similar masses are 
compared. Stassun et al. (2004) also report a similar result for a sample of young solar-like 
PMS stars, particularly for those stars with M < 0.5 M0. Preibisch et al. (2005b) find 
significantly lower X-ray emission for accreting as opposed to non-accreting PMS stars (as 
selected on the basis of Ca II triplet emission equivalent width), for the restricted mass 
range 0.3-0.5 M0. They also find that the non-accretors show a better-defined correlation 
between Lx and Lboi than the accretors, that the median value of Lx/Lboi is nearly a factor 
of three lower for accretors, and that there is a weak anti-correlation of the fractional X- 
ray luminosity with accretion rate. In looking for correlations of X-ray emission with the 
presence/absence of accretion disks, but for which there is no direct evidence of accretion, 
Preibisch et al. (2005b) find, somewhat surprisingly, that there are significant differences in 
Lx/ Lboi for stars with and without near-IR I — K excesses, but no differences for stars with 
and without K — L excesses; K — Lis usually considered to be a more reliable disk indicator. 

In our particular case, we have three different disk and/or accretion indicators (/ — 
H — K, and U — V excesses) for both clusters, and a fourth indicator (Ha equivalent widths) 
for NGC 2264. We note that the Ha measurements were only obtained for stars with 
classification spectra in NGC 2264. Most of these spectra were obtained as part of the 
studies by RebuU et al. (2002) and Makidon et al. (2004), where the goal was to classify 
stars with known periods. Therefore, the sample of stars with Ha is likely to be biased 
toward stars that are brighter in X-rays (see Figure 6). 

The numbers of stars available from the best possible sample are summarized in Table 3. 
While I — K and H — K excesses indicate the presence of either passive or actively accreting 
circumstellar disks, U — V excesses and Ha arise from active accretion; U — V excesses have 
already been translated directly to mass accretion rates for many of our targets (RebuU et al. 
2000, RebuU et al. 2002). 

Based on previously published results, we might expect that any correlations between 
total and fractional X-ray luminosity and stellar properties are likely to be weak (and in fact 
they are). We show plots for two of our available diagnostics of disks and/or accretion in 
Figures 10 and 11. Note that, for the near-IR disk indicators, redder colors suggest disks, 
whereas for U — V, bluer colors suggest accretion. The placement of the disk cutoffs is 
discussed in RebuU et al. (2002, 2000), along with comparison of the disk indicators for the 
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same stars. 

There is a very large amount of scatter in these two figures, plus their companions 
{H — K and Ha) which are not shown, although the scatter is somewhat smaller in the plots 
with log Lx/Lboi than in the plots against log Lx. For this reason, we have chosen to show 
only the box plot representations of these data. 

No statistically significant or consistent trends are seen in Fig. 10 or 11 (or in similar 
plots with the other available disk indicators). In stars without evidence for disks, we find 
examples of stars in which Lx is high and log Lx/Lboi is close to —3, but comparable large 
values are not seen in stars with disks or active accretion. The upper bound on both the total 
and fractional X-ray luminosity appears to decrease with increasing I — K excess for NGC 
2264 (but not the FF) and in increasing U -V excess for the FF (but not NGC 2264). The 
X-ray fluxes do not appear to depend on H — K excess for either cluster or in Ha for NGC 
2264. These confusing, inconsistent results suggest that if there is any trend, it is subtle, 
and, despite the numbers of points we have (see Table 3), we require many more data points 
to securely detect it. 

Figure 12 looks for correlations explicitly with log M as derived from the U — V excess as 
in RebuU et al. (2000, 2002). The data in the Orion FF are too sparse to reach any conclusion. 
For NGC 2264, however, there may be evidence of a trend: the upper envelope of L^/ Lboi is 
near 10~^ for stars with accretion rates < 10~^ M yr~^, but progressively decreases by nearly 
a factor of 10 for stars with higher accretion rates - a result in concordance with previously 
reported anti-correlations between X-ray luminosity and accretion for solar-like PMS stars. 

It is perhaps surprising that there are no clear correlations of X-ray emission with either 
the presence/absence of accretion disks, or the magnitude of accretion rate apparent from 



Table 3. Summary of numbers of stars available with and without disks 





Orion FF 


NGC 2264 


disk indicator 


non-disk disk 


non-disk 


disk 




detections (limits) detections (limits) 


detections (limits) 


detections (limits) 


I — K excess 


129 (38) 109 (37) 


226 (42) 


63 (21) 


H — K excess 


171 (56) 64 (15) 


236 (56) 


52 (15) 


U — V excess 


28 (4) 25 (3) 


87 (14) 


151 (48) 


Ha eqw 




77 (8) 


46 (12) 
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Fig. 12. — Log Lx and log Lx/Lboi vs. M (derived using DAM masses) for the best samples 
available for each cluster (see §2.3). There is a trend in L^, which is a rediscovery of the 
mass-dependent effects of found in Fig 4 and of M found in Rebull et al. (2000). There 
is no trend in Lx/I/bob but the dashed lines trace a trend with the highest M values. 
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our data, particularly given the plausibility of proposed X-ray production and absorption 
processes. Perhaps the intrinsic range in X-ray emission properties is so large that trends 
are difficult to discern, even with samples as large as several hundred stars. Alternatively, it 
may be that multiple mechanisms are at work simultaneously, and sorting among them may 
require a somewhat finer characterization of X-ray properties than Lx - the single property 
thus far used in searching for accretion-related correlations. 

We anticipate that the COUP data, in combination with the Orion Treasury program 
(Robberto et al. 2005) and the Spitzcr data in Orion (Megeath et al. 2006) will provide 
considerably more insight into the physical processes accompanying disk accretion. In par- 
ticular, the Treasury Program will provide robust estimates of ultraviolet excess emission for 
a large sample, which in combination with extant optical spectroscopic studies yielding R 
and Teff estimates, should provide accretion rate estimates for a large (perhaps 1000 stars) 
sample of accreting PMS stars. The deep COUP exposures of the ONC (20 x the length 
of the Chandra observations reported here) should provide well-characterized X-ray spectral 
energy distributions from 0.5 to 5 keV - sine qua non for sorting the relative contributions 
of soft X-ray emission that may arise from accretion, harder coronal X-ray emission and the 
effects of absorption arising in mass outflows. 



7. Conclusions 

We have studied the X-ray properties of stars in the Orion FF and NGC 2264. With 
ages of '^l Myr and ~3 Myr respectively, both regions contain stars that are older than 
most stars in the ONC. Rotation periods are available for a relatively large fraction of the 
PMS population in these regions. Furthermore, NCC 2264 contains a significant number of 
stars that have made the transition to radiative tracks. These samples thus enable study of 
evolution- and rotation-driven changes in X-ray properties over the age range 1-5 Myr. 

The data for the Orion FF and NCC 2264, when combined with the much more extensive 
data from the COUP survey of the ONC, establish a number of clear patterns in the behavior 
of X-ray emission in PMS stars. In all three regions, log for stars on convective tracks 
correlates positively with both log Lboi and also with stellar mass, not surprisingly since 
log Lboi is strongly correlated with mass. As Preibisch et al. (2005b) point out, a strong 
correlation between X-ray emission and stellar properties argues that the X-ray emission 
probably originates primarily in compact (coronal?) structures associated with the star. We 
confirm as well that the level of X-ray emission is directly related to stellar structure in that 
the level of emission drops on average by ~10 when PMS stars with masses in the range ~ 
1-2 Mq make the transition from fully convective to radiative tracks, where stars of about 
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solar mass have radiative cores and convective envelopes. 

The fractional X-ray luminosity, log L^/ Ly^oh is not strongly correlated with mass for 
stars on convective tracks; the scatter in the relationships of both log Ly_ and log Lx/Lboi with 
mass is large. The physical reasons for this large scatter remain unclear. Age is apparently 
not a large contributor; while there is evidence the ONC has nearly twice as many stars with 
extremely high X-ray emission (log Lx/I/boi ^ —2.9) as the Orion FF and NGC 2264, we 
see little change in median values of either log or log Lx/Lboi during the first ~5 Myr 
of evolution down convective tracks. In sharp contrast to main sequence stars, we find no 
correlation between log and log L^/ Lboi with P oi v sin i for stars on convective tracks, 
although more vsmi data arc needed in NGC 2264. Preibisch et al. (2005) have argued 
that X-ray variability at typical levels is likely to be too small to account for the scatter in 

-^x/ -^bol- 

In both the Orion FF and NGC 2264 we find, in agreement with earlier studies of the 
ONC, that stars with larger values of log Lx/Lboi are more hkely to exhibit spot- modulated 
periods. Approximately 3/4 of the stars with log Lx/I/boi> ~3 are periodic, while while 
only about 1/3 of the stars with log Lx/I/boi< —4 have periods, and this difference does not 
appear to be the result of any biases in the samples. It may well be that, as Stassun et al. 
(2004) have suggested, more active stars have larger and/or more organized magnetic spot 
coverage and that periodic variations are more easily measured in such circumstances. 

We have searched for a correlation between log Lx and log Lx/Lboi and 1) near-IR excess, 
which is a circumstellar disk indicator; 2) UV excess and Ha, both of which are indicators of 
accretion; and 3) the mass accretion rate, as inferred from the UV excess. There is substantial 
scatter in all of these relationships, and no clear trends emerge. Discovery of the relationship, 
if any, between accretion and X-ray emission awaits the combination of COUP results, which 
provides measures of X-ray luminosity and X-ray spectral energy distributions, with the 
corresponding results from the ongoing HST Orion Treasury program, which promises to 
provide accurate estimates of accretion rates for large samples of PMS stars, and the Spitzer 
data in Orion, which will provide more accurate mid-IR disk indicators. 
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